CORRELATIONS FOR THE CIRCULAR DYSON BROWNIAN 
MOTION MODEL WITH POISSON INITIAL CONDITIONS 



P.J. ForresteiQ and T. Nagac0 



The circular Dyson Brownian motion model refers to the stochastic dynamics of the 
log-gas on a circle. It also specifies the eigenvalues of certain parameter-dependent 
ensembles of unitary random matrices. This model is considered with the initial con- 
dition that the particles are non-interacting (Poisson statistics). Jack polynomial 
theory is used to derive a simple exact expression for the density-density correlation 
with the position of one particle specified in the initial state, and the position of 
one particle specified at time r, valid for all (3 > 0. The same correlation with 
two particles specified in the initial state is also derived exactly, and some special 
cases of the theoretical correlations are illustrated by comparison with the empirical 
correlations calculated from the eigenvalues of certain parameter-dependent Gaus- 
sian random matrices. Application to fluctuation formulas for time displaced linear 
statistics in made. 



1 Introduction 

The Dyson Brownian motion model refers to the overdamped Brownian dynamics of the 
one-dimensional log-gas. This dynamics is specified by the Fokker-Planck equation 



Cp where -C = V ^— ^ h 

oxj \dxj dxj^ 



dr 

with the particular potential 



(1.1) 



or its periodic version 



1 ^ 

W = -Y.^]- E log\xk-xjl (1.2) 



W = - Yl logle^™'--/^ -e^^^^^/^l. (1.3) 
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The formulation of this model has its origin in the theory of parameter dependent Gaussian 
random matrices. We recall that a random matrix H is of this type if its joint distribution of 
elements is proportional to 

exp ( - pTr{H - e-^H^^'> f/{2{l - e-^"))) , (1.4) 

where /3 = 1,2 or 4 according to H being symmetric, Hermitian or self dual real quaternion 
respectively. This means that each part (real or imaginary) of each independent element of H is 
chosen with a particular Gaussian distribution. The matrix H^^^ is some prescribed (random) 
matrix. Note that as r ^ 0, H = while as r — > oo dTD reduces to exp{-PTvH^ /2). 

With P = 1, 2 and 4 this latter distribution specifies the Gaussian Orthogonal, Unitary and 
Symplectic Ensembles (GOE, GUE, GSE) respectively. 

Let p = p(xi, . . . ,XAr;r) denote the eigenvalue probability density function (p.d.f.) corre- 
sponding to ([L.4|). By using second order perturbation theory, Dyson ^ showed that p satisfies 



the Fokker-Planck equation (LI) with W given by (L2). This equation must be solved subject 
to the initial condition that p corresponds to the eigenvalue p.d.f. of H at t = 0. Dyson also 
formulated a parameter dependent theory of random unitary matrices from the so-called circular 
ensembles. Here the situation is more abstract in the sense that no explicit construction of such 
unitary random matrices is known. Nonetheless, in the framework of the abstract formulation 
it was shown that the corresponding eigenvalue p.d.f. (eigenvalues ^"^^^il^ ^ < Xj < L) satisfies 
the Fokker-Planck equation (0) with W given by (pT^). It turns out that upon introducing the 
scaling 

X ^ irpx/Vm, T ^ {■KpfT/{2N) (1.5) 



into the Fokker-Planck equation with W given by (|l.2| ), and taking the N ^ oo limit, the results 
obtained for the correlation functions and other observable quantities are identical to those 
obtained with the choice of W ( |1.3D . In fact for correlations over one or two eigenvalue spacings 
in the bulk of the spectrum, the agreement already becomes apparant for matrix dimensions 
= 11 (see e.g. [§). 

Our specific interest in this paper is the exact calculation of certain correlation functions 
in the Fokker-Planck equation ( |1.1| ) with W given by ( |1.3| ) and subject to the initial condition 
p{xi, . . . ,xn;0) = pv. This initial condition corresponds to non- interacting particles; in the 
applied random matrix literature it is referred to as Poisson initial conditions, which is consistent 
since a gas of non-interacting particles exhibits Poisson statistics. Thus for the random matrix 
couplings /3 = 1,2 and 4 and N large, the Fokker-Planck equation with this initial condition 
describes the transition between an eigenvalue p.d.f. with Poisson statistics to the eigenvalue 
p.d.f. of the appropriate Gaussian ensemble. Such a transition is relevant to the description 
of the statistical properties of the eigenvalue spectrum in a quantum system which is initially 
integrable, but becomes chaotic as a parameter (usually identified with r^/^; see e.g. Q) is varied. 
Because of this application, this problem has received a lot of recent attention |^, |6|, |^, ^ . 

In the next section we will specify the particular correlation functions to be calculated 
(density-density correlation with the position of n particles specified in the initial state, and the 
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position of one particle specified at 'time' r), and then proceed to calculate this correlation in 
the case n = 1 for general (3. Comparison with the empirical evaluation of this correlation for 
(3 = 1 obtained from numerically generated parameter dependent random matrices is also made. 
In Section 3 we consider the same correlation with n = 2, and provide its exact value for (3 = 2 
and 4, as well as all < /3 < 2. A discussion of these results, including asymptotic properties 
and their relationship to fluctuation formulas for linear statistics, is given in Section 4. 



2 Density-density correlation 
2.1 Formalism 

In general the correlation functions for the Brownian motion described by ( |1.1| ) can be specified 
in terms of the Green function G{x^i \ . . . , x^^; xi, . . . , xat; r), which is by definition the solution 



of (LI) subject to the initial condition 



TV 

p(xi,...,x„;0) = \[5{xj-xf). (2.1) 

i=i 

For a general initial condition with p.d.f. f{x^i \ ■ ■ ■ (assumed symmetric) the particular 

density-density correlation ^■^{x^i \ . . . , Xn^; x; r) is given in terms of G by 

pf„,i) (xf ) , . . . , X W ; x; r) = iV(iV _ 1) . . . (AT _ ^ + 1) dxf^, ■ ■ ■ dx^J^ f{xf^ , . . . , x^J^ ) 

X / dx'i^ ■■■ dx^^ {^^{x - ^f^ijGixfK. . . ,x^^;xi, ...,xn;t) 
Ji Ji 

-p(„)(xS°\...,xW;0)p(i)(x;r). (2.2) 

Note that ^^(x^''\ . . . , x^?'*; x; r)/p(„)(xj^''\ . . . , Xn^"*; 0) -|- p(i)(x;r) represents the density at 
position X after time r, given that there are particles at x^*^"*, . . . ,Xn^ initially. 



To obtain a useful expression for the Green function, one uses the general fact |10| that after 
conjugation with e~^^ the Fokker-Planck operator C transforms into a Hermitian operator. In 
fact (see e.g. |ll|) for the choices of W (|1.2|) and (pi) one has 



where H is the Schrodinger operator for a quantum mechanical system with one and two body 
interactions only, and Eq is the corresponding ground state energy. Explicitly, for W given by 



which is an example of the Calogero-Sutherland quantum many body system (quantum particles 
with pair interaction). Now in general for the imaginary time Schrodinger equation 

^ = -;g(^-^o)V (2.5) 
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in which H has a complete set of orthogonal eigenfunctions {iPk{x)}k (here x := (xi, . . . ,xj\f) and 
K represents an n-tuple of ordered integers ki > k.2 > ■ ■ ■ > k^) with corresponding eigenvalues 
{Efi}f^, the method of separation of variables gives that the Green function solution is 

In (P) 

|2 



(V'kIV'k) := dxi--- dxN\^IJn{x)r, (2.7) 
and * denotes complex conjugation. 

Since, according to ( |2.3D , the Green function solution G{x^^^ ; x;t) of the Fokker-Planck 
equation ( |L1] ) is related to the Green function solution G^^\x^^^;x;t) of ( ^^ ) by 

p-/3W^(xi,...,xjv)/2 

G(x<»);x;.)= ^^^^^ gW(...);.;,). (2.8) 



G{x^^^ ; x;t) is determined to the extent that the quantities in (2.6) are known. Independent of 



such explicit knowledge, substituting (^^) in ( |2.8| ) and substituting the result in ( |2.2| ) gives 

p^^,^ (xf ) , . . . , X W ; r) = iV(iV - 1) . . . (iV - n + 1) (2.9) 
, V -4.(xf \ . . . , xg'))(e-^^/^| Ef=i ^(x - x,)|V^.) 



where 



yl.(xf\...,xW) = y^dxi°j,...y^dxSV(xS°\...,xW)V^.(x(°))e/^^/2 (2.11) 

TV 

(e-^^/2|^5(x-x,-)|^.) = iVy^(ix2---y^dx^e-^^(^')/2^,(x). (2.12) 



2.2 Explicit formulas 

For the Schrodinger operator (p.4|), it is known |12, 13 1 that the (unnormalized) eigenfunctions 
can be written in the form 

^,{x) = e-^^/2p(2//3)(^)^ ^ _ g2Wi, (2.13) 

where k = (ki, ... , kjv) forms a partition and thus has all parts non-negative and ordered as 
specified below and P^^^^\z) is a particular polynomial known as the Jack polynomial. 



Each Jack polynomial has the special structure 

(2.14) 



tJ.<K 



where refers to the monomial symmetric function (in the variables zi, . . . , zn) corresponding 
to the partition k, the b^fj, are coefficients, while n < n refers to the dominance ordering of 
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partitions. For k possessing negative parts, ( |2.14| ) has no immediate meaning. In such cases the 
Jack polynomials are defined by the relation 

N 

P!.%'\z) = {U-r)Pi'^'Hz), (2.15) 
1=1 

where k is a partition while c := (c, . . . , c), c G 2^. Also, for general k, the energy eigenvalues 
are given by 

E,-Eo = {-) (5:4 + |^K,(iV-2j + l)). (2.16) 



Although (|2.14| ) is a very special property associated with (p.4|), the true utility of Jack 
polynomial theory in regard to computing correlation functions lies with our knowledge of the 
integrals ( |2.7|) and ( |2.12| ) (see e.g. [^] and references therein). First consider the normalization 
integral ( |2.7D . In the case of the ground state (k = 0) we have 



(V'olV'o) := [ dxi--- [ dxN T] \e 

Jo Jo ,^jtt^AT 



2mXj/L g27riZfc/L 1/3 



^<j<k<N 

_ ^ n(iN/2 + 1) 

" (r(/3/2 + i))^- ^'■"> 

This can be conveniently factored out of the normalization integral for general k. In addition, 
the general k case involves the generalized factorial 



N r(n-f(j-i) + K ,) 
,A r(^-f(i-i)) 



^ -^-7^;-;^— ^^ (2.18) 

the magnitude of the partition, 



N 

\K\:=Y.Kj, (2.19) 

i=i 

the Jack polynomial which each variable set equal to unity pj^^^\l^), and the quantity 

<■= n (^•-^ + |('^^-i + l)). (2.20) 

In ( ^) it is necessary that k be a partition The product is over all squares in the diagram of 
K, and k! refers to the partition conjugate to k. In terms of this notation, for k a partition we 
have 

mi^.) = mi^o r' ^ ^ IJ - (2.21) 

[|(iV_l) + l](2//3) 

The normalization of the states in which k possesses negative parts and so does not form a 
partition are reducible to the latter case by the simple formula 



N N 

n^rjv'«|(n^rjv'K) = ii^M- (2.22) 

1=1 1=1 
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We turn now to the integral ( 2.12| ). This quantity has the evaluation 



N 



1^1(^1-1)! 
L 



N 



na(-fo--i))^^. 

[f(iv-i)+ir^) ' 



(2.23) 



where ^(k) denotes the length of the partition k (i.e. the number of non-zero parts) and 

(a)„ := a(a + !)••• (a + n- 1). (2.24) 

For the n-tuples of non-positive integers 

- K := (-KAf, -KAT-i, . . . , -Ki) (2.25) 

we have the simple formula 

N N ^ 

(V^o|E'^(^-^i)|V'-K) = (V'o|E'^(^-^j-)|V'k)*, (2.26) 

j=l 3=1 

and thus (|2.23| ) gives the evaluation. Furthermore, for /3 rational, we have the result that 
( 2.12| ) is non-zero only if all parts of k are non-negative or all parts are non-positive. Hence, if 
we restrict attention to /3 rational (which we will do henceforth), the results ( 2.23| ) and ( p. 26 ) 
suffice. In this regard we note also that it is a simple result that 

{^--M-r) = ii^M, (2.27) 



which when combined with ( 2.26 ) implies that we can restrict the summation in (|2.9| ) to parti- 
tions provided twice the real part is taken. 

Finally we come to consider the integral defined by ( 2.11| ). For the case under investiga- 
tion, this reads 

which is proportional to the term independent of z^^-^, . . . , in the power series expansion of 
Since we have now reduced the cases to be considered down to those for which k is 
a partition, this term can be determined by setting z^l^^i = ■ ■ ■ 



,(0) 



and we therefore have 



(0). 



L 



(2.28) 



,zf,0,...,0)) 

otherwise, 

where the first expression in the second equality follows from the fact that the coefficients ft^/x 
in ( p.l4| ) are independent of N. 

Substituting ( |2.2l|) , (|2.23|) and (|2.28| ) into and taking twice the real part we obtain the 
formula 



N{N -1)---{N -n + l) 



(0). 



(2.29) 



partitions k, 
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where 



L <(/3/2)l"l 2 



vahd for /3 rational at least. In fact this result must be valid for all /3 > 0, by continuity of pj^ ^-^ 
in p. 



2.3 The case n = l 

In the case n = 1 the only partitions contributing to ( |2.29| ) are (ki, 0, . . . , 0) and we have 

Pi'^'Hz?') = E.-Eo = {^y{.l + - 1).,), = 1. (2.31) 



Thus ( |2.29| ) reduces to the remarkably simple formula 



T .^J0)__, _ 2(A^-1)P ^ e-(2-A)'('^?+(^/2)(^-i)'^^W^cos27rKi(x-xS°V^ (2-32) 

K1=0 



where p := A^/L. In the thermodynamic limit N,L ^ oo, p fixed, ( p. 321) being a Riemann sum 
tends to the definite integral 

pfi,i)(xS°);x;T) =2/,2 J^-i^-P^'i^'+W^^^^-Z^ cos27Tspix - x^^"^) ds. (2.33) 



As mentioned in the Introduction, the correlation function given analytically by (|2.32| ) with 
P = 1,2 or 4 can be calculated empirically using parameter-dependent random matrices. For 
definiteness consider the case /? = 1. We construct initially a diagonal random matrix H^^^ = 
diag(/i^'J\ . . . , /i^]y) with each diagonal entry chosen from a Gaussian distribution with mean 
zero and standard deviation A^(27r)^^/^. This standard deviation is chosen so that the eigenvalue 
density, which equals times the Gaussian distribution for a diagonal element, is unity at the 



origin. Then, following the prescription ( |1.4| ) with /? = 1 and the scaling (1^), we construct a 



real symmetric random matrix in which the diagonal entries have mean e.~^h^jj and standard 
deviation (\/2A^/vr)(l — e~^*)^/^, t := 7t'^t/(2N), while the independent off diagonal elements 
(upper triangular elements say) have mean zero and standard deviation equal to 1/ \/2 of that 
of the diagonal elements. Note in particular that the factor of y/2N /vr in the standard deviation 
is chosen so that at the centre of the spectrum the theoretical density remains equal to unity. 

We choose a specific value of and t, and numerically generate many such random matrices 
with a fixed initial diagonal entry h^^^ = 0, together with their eigenvalues. The corresponding 
empirical eigenvalue density in the neighbourhood of the origin can be computed and compared 
directly against the theoretical result ( p.32| ) with p = 1, /3 = 1 and the equilibrium density p = 1 
added. The results of such a calculation are given in Figure |l[ We see that the empirical and 
theoretical results agree to statistical accuracy. 
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Figure 1: Comparison between the empirical density for 2,500 Gaussian real symmetric 
parameter-dependent random matrices of dimension 13 x 13 with t = .025 and a zero 
eigenvalue initially, and the theoretical value computed from ( p.32| ) with /3 = 1. 



3 The correlation p^^^)(xi°\ ^2°''; r) 

In this section we will consider the correlation ( p.29 ) with n = 2. Explicit formulas are possible 
in this case because of the formula [^ ] 

ptJ/lM'\4'') = ^(4°^4°V--^-)/^^^t (^(4°^ + )-/^), (3.1) 

where on the r.h.s. Pn°''^\x) denotes the Jacobi polynomial, 7 := (/3 — l)/2 and 

a.:= r + 7) (- + ^; + l)- 2-". (3.2) 
V n / (7 + l)n 



Of particular relevance is the asymptotic form of ( p. ID in the limit ki,K2,L — > 00, p fixed. 
The leading behaviour of the term (|3.2| ) is determined using Stirling's formula, while the Jack 
polynomial is estimated by noting that 

+ )(.S°).f )-/2 = eosvr(xf) - (3.3) 

and making use of the asymptotic formula 

Pj{^'^\cos9)^Qyj,{ne), (3.4) 

where J-y denotes the Bessel function. Thus we have 



where si := Ki/N, S2 ■= K2/N . 



8 



Due to the factor {-(3/2) 

K2 111 

(l2^ ), the cases for which /3 is even require special consider- 
ation since then K2 is restricted to the range to (3/2 — 1 for a non-zero contribution. In fact, 
due to cancehation effects within this range, each even value of (3 must be considered separately. 
We will consider in detail the cases (3 = 2 and 4. 



3.1 (3 = 2 



In this case we see from ( p. 30 ) and ( 2.20 ) that with K2 = 0, = \/L while for K2 = 1, 
Un = — 1/L. This fact, together with (|3.5| ) and ( 2.16| ) shows 



1/2 



K2=0 



xe 



^{2^p)2{.2+si)r/2f_^ 



1/2 



-M ._ ^(0) ^ ^(0) ._ ^(0) 



1/2 



Ji/2(^^ip(4°^-xf ))}, (3.6) 



where x\ := x\ — x, X2 := x\ — x. Making use of the formula 

/ 2 n1/2 

•^1/2 (^) = ( — sin 2;, (3.7) 
summing over k\ and substituting the result in ( p. 29 ) gives that in the thermodynamic limit 

pl2,x)^xf\xf-x-T) = 2p3 ds cos7rsp(xS°) +xf )e-(2-'')'(^'+^)-/2 

(0) ^;^(0)^ -^^^/^.W ^^(0)^ 



t (0) (0)x 

—T^^ ^y— ) [^P{x\ + x\ ')sin^sp{x\ + x\ 

■np{x\ ' -x\') 

(0) , ™(0)^^p-{2^p)'{^'+^)r/2 



As an analytic check on (|3.8|), consider the limit ^2^^^ 00. Rewriting the trigonometric 
products as sums shows that 



+2p-' / ds 



+{2-Kp) —cos'Ksp{x\ + X2 ))e. 



(3. 



T /^(o)_.. 



(0) 
— >oo 



(3.9) 



as expected 
A case ( 

state. Making use of the formula 



A case of special interest is when x^^^ — x'2' ^ so that two particles coincide in the initial 



(3.10) 



in (p. ID and then proceeding as in the derivation of (3.5) we have that in the finite system 



T I (0) (0) N 
P(2, 1)1^1 '^1 ^X^T) 



2{N - l)(^^-_2)j^^ g-27rmip{s;("^-:r)/Arg-(27rp/Ar)2{K2 + (Af-l)Ki)r/2 



L3 



(( 



x((ki + 1) - Kie 



Kl=l 

2TTip{xf'^ -x)/N -{2-Kp/Nf{N-2)T/2 



(3.11) 
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1.8- 




Figure 2: Comparison between the empirical density for 2,000 Gaussian Hermitian 
parameter-dependent random matrices of dimension 13 x 13 with t = .05 and a doubly 
degenerate zero eigenvalue initially, and the theoretical value computed from ( p.ll|) . 

It is easy to see that in the thermodynamic hmit ( |3.11| ) agrees with ( |3.8| ) with n , thus 

providing another check on the latter resuh. We can also compare the correlations obtained 
from the analytic formula ( |3.11| ) with those obtained empirically from parameter-dependent 
Hermitian random matrices with the initial condition of a doubly degenerate zero eigenvalue. 
The parameter-dependent matrices are constructed in an analogous way to that described at 
the end of the previous section. Thus the initial matrix H^^^ is again diagonal and real, with the 
elements hj^j, j > 2, chosen with the same Gaussian distribution as before, while the elements 
/ill and /i22 are set equal to zero. The matrix H is chosen so that the diagonal entries have mean 
e^^hfj and standard deviation (\/]V /7r)(l — e~^*)^/^, t := 7r^r/(2A^), while the independent off 
diagonal elements (real and imaginary parts of the upper triangular elements) have mean zero 
and standard deviation equal to l/\/2 of that of the diagonal elements. Again one feature of 
this prescription is that the theoretical density at the origin is unity. 

In Figure |2| we give a comparative plot of the empirical density obtained from generating the 
eigenvalues of such matrices for a certain choice of N and t with the corresponding theoretical 
formula ( 3.11J ). 



3.2 [3 = 4 



Here the formulas (|2.30|) and ( ^.20|) give 



K2=0 



I' 



K2 = l 



2 / 



1 



Ki + 2 



K2=2 



Kl + 1 



(3.12) 



These values suggest we write the corresponding integrand for the sum over ki in the form of a 
difference approximation to the second derivative. Doing this, and making use of ( |3.5D we find 
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that in the thermodynamic Umit 

x(£/(^) - (-/2)(2-p)V(^) + 2|(^)) (3.13) 

where 

f{s) := s^H ^_(o) - #).)e-(^i"+4"')^e(-'')^-/2. (3.14) 

Regarding checks on this formula, from the large-x expansion J3/2ix) ~ — (2/(7rx))^''^ cos x 
it is straightforward to show that the hmiting behaviour (|3.9| ) is exhibited. Furthermore, use of 
( |3.iq ) as well as ( p. 121 ) in ( |2.29D gives that in the finite system with Xj^^ — 



2 ; 



3 

Kl = l 



-2(ki + l)2e-(27rp/iV)2(l+2(Af-3))T/4g27ripx('"/A^ + (^^ _ ^_^^-i27Tp/Nf {4+4(N-3))t/4 ^4nipx^°'> /Ny 



(3.15) 



In the thermodynamic limit this agrees with (|1|) in the case x^i^ = x^^ (the formula Ju{x) ~ 
x'^/(2'^r(l + I/)) as X —> is required in the checking). 

3.3 General < /3 < 2 

In the cases (3 = 2 and /3 = 4 we have seen that cancellation takes place within the summand of 
( 2.29| ) due to the factor (— /3/2)k2 having varying sign. This is again true for general < /3 < 2. 



Thus for = 0, (-/3/2)«2 = 1, while for K2 > 1, (-/3/2)«2 = -(/?/2)r(K2 -/3/2)/r(l-/3/2) < 0. 
However in the case of general < (3 < 2 the factor {—(3/2,)^^ is non-zero for all K2 > 0, so the 
method of grouping together terms for cancellation used in the cases /? = 2 and /? = 4 is no 
longer applicable. 

To gain some insight into the necessary grouping, let us consider the behaviour of as 
specified by ( 2.3C| ) with £{k) = 2, for large ki, K2- Use of ( |2.20| ) and Stirling's formula shows 

iP/2) (^Ji + K2)(k1 - 1)! r(K2 - P/2) TiP/2 + KI-K2 + 1) 



L K2\{Ki-K2y- r(l-/3/2) r(/3/2 + Ki + l) 
{(3/2) {K^ - K2f/^ 



LT 



(3.16) 



{l-P/2) ^l+/3/2^1+/3/2- 

Writing this in terms of continuous variables si := ki/N, S2 '■= K2/N we see that ( p.l6D is 
non-integrable at si = and ,82 = 0. It is this singularity that we must cancel out using the 
K2 = term. 

We proceed as follows. Define 

_ (ki + K2){ki - 1)! T{P/2 + KI-K2 + 1) 

"'^''"'''^ L{ki-K2)\ r(/3/2 + Ki + l) 
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so that 



and 



iP/2)T{K2-P/2) 
r(l-/3/2)K2! 



(3.18) 



(0) jo)__,_iV(iV-l) 
L2 



;2Re( ^ ^;(^^,o)^(«l,o)(4°^4°^) + fi'K2^(«i,K2)^(«i,«2)(4°^ 4°^)) > (3-19) 



Kl = l 



Kl>ft:2>l 



where 



Now write 



B. 



fi'K2 — G^a+i — Gk2, Go :— 0, Gi :— 1. 



From the general formula 



N 

(afc(^A;+i — fefe) + fefc+i(afc+i — fflfc)) = OAT+i^iv+i — ao^O) 



(3.20) 
(3.21) 

(3.22) 



fc=0 



we see that ( 3.19 ) can be rewritten 



T / (0) (0) X 
P(2,1)V'^1 ' -^2 T^j'^) 



N{N - 1) 

l2 



x2Re( J2 G'«,2(^(«lA2)^(«lA2)(4°^4°^) -^{«iA2-i)-S{«i,«2-i)(4' 



(0) (0)- 



; ^2 



'il>K2>l 



(0) JO) 



Kl = l 



(3.23) 



The utility of ( 3.23| ) is that the summand is well behaved in the large ki, K2 limit. To see 
this, note from the definition ( |3.18| ) that for large K2, 

(/3/2) 1 



9k2 



r(i-/?/2) 



and so ( |3.21 ) is asymptotically satisfied with 



Gk2 ~ 



r(i - p/2) ■ 



(3.24) 



(3.25) 



This has an integrable singularity at the origin, as distinct from the non-integrable singularity 
in ( |3.16| ). Making use of the asymptotic formulas ( 3.25 ) and (|3.5| ) and noting from ( 3.17 ) that 



K ~ ^(k1 + K2)|ki - K2|'^/^^ 



we see from ( p. 23 ) that in the thermodynamic limit 



T I 



(0) (0) X 
1 I •''1 ' J 



1+/3/2 ■ 



(/3-l)/2 ^oo /-oo 

dsi I ds2 



(/3/2)2r(-/?/2) V27rp(xf) - 4°V h ' Jo 



(3.26) 



1 \/3/2 



SlS2 



X ((Sl + S2)\S, - «2|(^+^)/V(,_i)/2(Ni - S2\np{xf^ - #)) 

X COS (vrp(xf) + 4°))(si + ,2))e-(2-^)'(^?+^^+(^/^)(^^+^^))-/^) , (3.27) 



where an integration by parts in si has been carried out. 
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4 Discussion 



4.1 Asymptotics 

In the guise of the Calogero-Sutherland model, the density-density correlation pj-^ -^-^{x^^ ; x;t) 
for the Dyson Brownian motion model in the case that the initial state equals the equilibrium 
state (in.= eq.) is known for all rational (3 (see or [l^ for the explicit connection between 
the correlations of the two models). For P/2 = p/q (p and q relatively prime) its value in the 
thermodynamic limit is given in terms of a (p + (7)-dimensional Dotsenko-Fateev-type integral. 
This contrasts with the simple one-dimensional integral (2.33) for the same correlation with 
Poisson initial conditions. 

Although the expression for pj-^ -^^{x^^ ; x;t) in the case in.= eq. is complicated, its non- 
oscillatory leading order large \xj^^ — x\ and/or r expansion is very simple, being given by 



4p2 



pJ...{x^!'^;x;t) ~ ^Re( ^). (4.1) 

'^(i.i)^ 1 ^n.=eq. (3 \ i27r p)^ + 27ri pxf^ )^ ' 



Fixing r, ( [4.1| ) shows that the leading non-oscillatory portion of -|^^(x^'^^; x; r) falls off as 
— l/{/3{7rxi'^)'^), while with xf'^ fixed the decay l/(/3(7r yo)^r)^ is exhibited. Note that the decay 
in xf*^ is independent of the density. 



Let us compare (4J) with the corresponding asymptotic expansion in the case of Poisson 
initial conditions (in.= Poi.). By rewriting the cosine term as a complex exponential, linearizing 
the exponent about s = 0, and changing variables we see from ( |2^ that 



^a..,(^r';:.r),^^^^^^VRe^ ^^^^^^^;^^^^^,„, J ,4,2) 

Notice that there is no /3 dependence in this expression. For xf'^ fixed the decay here is given 
by l/(7r^T), independent of p, while for r fixed the decay is p^T / {x^^^ . These behaviours have 
distinct features from those noted above for ( [4.11 ). 

The leading order large xf'^ and/or r asymptotic expansion of ( |3.8| ) {(3 = 2 result) with 
x^^^ = x^2^ can readily be computed. We find the same behaviour as in ([4.2|) , except that the 
prefactor is replaced by Ap^. For general < (3 < 2, setting x^'^^ = in (|]23), then 
repeating the analysis which led to ([4.1| ), we find 

Ph) ' ' ^' ^) ~ V^(/3)Re (- ^) (4.3) 



^{271 p)'^ - 2mpxY 



where 

m 



7rV2 1 



2/3+i(/3/2)2r(-/3/2)r(/3/2 + l/2) 

X dsj ds2( ) j^^({si + S2)\s,-S2fe-^'^+''^). (4.4) 

70 Jo ^SiS2^ dsids2^ ' 
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The integral in ( [4. 4]) is evaluated in the Appendix. Subtituting its value gives A{P) = 1, so the 
asymptotic behaviour found at /3 = 2 persists independent of /3 (for /3 < 2 at least) analogous 
to O. 



4.2 Fluctuation formulas 

Let us now turn attention to the application of the density-density correlation pj^ {x^^^ ; r) in 
the study of fluctuation formulas. In general [^] the time displaced covariance of two linear 
statistics 

N N 

Ar = J2aix,iT)), Br = Y^h{x,{T)) (4.5) 
is given in terms of the Fourier transform 

S{k-T) := r pfi^i)(xS°^;x;r)e^^i"'^dxf) (4.6) 

(assuming a fluid state so pj^ ^-^{x^^^ ; x]t) is a function of x^^^ := x^^^ — x) according to the 
formula 

/oo 
a{k)b{-k)S{k;T)dk. (4.7) 

Now, from ( p.33 ), for the Dyson Brownian motion model with Poisson initial conditions 

5(A;;r) =pe-'='"/^-"''l'=l", (4.8) 



so ( |4.5| ) is known explicitly. 

Rigorous studies of the infinite Dyson Brownian motion model for j3 = 2 and with in.= 
eq. , and of the same model on a circle ||2^ for general /? > have shown that after appropriate 
scaling the joint distribution of {Aq^Bj-) is a Gaussian. In the infinite system the covariance is 
given by ( [4.7D with S{k] r) replaced by its scaled form, while for the circle system, a discrete 
version of (^]^) applies. Explicitly, in the infinite system a small parameter e is introduced in 
(4.5) so that the linear statistics become 

N N 

A^ = Y, a{eXj{T)), Br = Y. HeXjir)). (4.9) 

Physically this means that the variation of a and b is macroscopic. Also, r is scaled by writing 

t = r/e. (4.10) 

Then it is proved in |19] that in the limit e — > and with /? = 2 the joint distribution of {Aq, Bt) 
is a Gaussian with covariance 



1 roa 

— / a{k)h{-k)\k\e-\^\^PUk. 
7rp J-oo 



14 



Also relevant to the present discussion is the corresponding result for the Brownian dynam- 



ics specified by (1.1) in which the equilibrium state is compressible (gas with a short range 
potential). Here one introduces the scaled linear statistics by 



N N 

Ar = ^ a{eXj{T)), Br = e'^^ ^ 6(ex,(r)) 



and scales r by writing 



t = r/e^ 



(4.11) 



(4.12) 



In this setting it is proved in |21| that the joint distribution of {Ao,Bt) is a Gaussian with 
covariance 



X / d{k)b{-k)e-^'p'/^^^'^ dk, 



where x denotes the compressibility. Note that as well as the different scaling for r, the linear 
statistics are suppressed by a factor e^/^ which is not required in ( |4.9|) for the Dyson model. 
This means that in the Dyson model the fluctuations are naturally suppressed by the long-range 
nature of the pair potential. 

In the present work we are considering the Dyson model with Poisson initial conditions. Thus 
initially the particles are non-interacting so the system is compressible. A simple calculation 
using ( [4.7p and ( |2.33| ) shows that if we scale the linear statistics according to (4.11) as in a 
compressible gas, but scale r according to ( fl.lO ) as in the Dyson model with in.= eq., then the 
covariance in the limit e — > becomes equal to 

l{k)b{-k)e-''p\''\^ dk. (4.13) 

We conjecture that in this limit the joint distribution of {Ao,Bt) is a Gaussian. 

Following [^] we can also consider the covariance for the system on a circle. For this purpose 
one first scales xf'^ and x by multiplying each by A^, and also scales r by making the replacement 
( 4.10| ). Then, according to ( 2.32| ), in the thermodynamic limit 

J2 oo 



2(7rp)2 \t 2TTip(x^°^ -x)k 



(4.14) 



where now x^i\x € [0, Ip]. In this setting, analogous to ( [4.11 ) we choose for the scaled linear 
statistics 

-j^ Af 1 ^ 

^- = 1^E«(^^)' Br = ^Y.^{x,). (4.15) 



Ari/2 . 

The formula for the covariance is then 



Cov(ylo,^r) 



1 



dxf'' a{xf'') I dxb{x)p^^i i^{xi" ,x), 



Ari/2 



T 



(4.16) 



where the factor 1/A^ results from the factors of 1/A^^/^ in (4.15), while the factor A^^ results 

^ A^x. Substituting (|1|) in (|4T^ ) gives that as 



from the change of scale x^ 

N ^ oo 



(0) 



Nx 



(0) 

1 1 



Cov(^o,-Br)= 



-2{7rp)2|Ki|t 



(4.17) 



K]^ — — OO 
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where 

/•Vp 

a^, := p / a(x)e2^*^^'^i dx 
Jo 

and similarly the meaning of ft-K^. Again we would expect the joint distribution of {Aq, Br) to 
be Gaussian. 
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Appendix 

Here we will evaluate the integral 

roo POO 1^2 



which includes the integral in (4.4) as a special case. Assuming temporarily that Re(;u) > — 1, 
integration by parts gives 

/•OO /"OO 

= (/i+l)2 / dsi dS2{siS2Y{si+S2)\si-S2\''e-'^'^+'^\ (A.2) 

JO JO 

Furthermore, the term (si + S2) can be obtained from e"^^*!"'"'^^) i^y partial differentiation with 
respect to z. Doing this, then scaling out the z dependence by changing variables and computing 
the derivative gives 

lip, z, (3) = {p + 1)2(2 + 2p + /3)z-(3+2m+/3)l(^, /3) (A.3) 

where 

/■OO /•OO 

Hp^P) := ds^J^ ds2{siS2r{si+S2)\si-S2fe-^'^+'^^ 

r(i + /.)r(i + /3)r(i + /i + /3/2) 

r(l + /3/2) ^ ■ ^ 

where the second equality in (A.l) follows because L{p, f3) is a two-dimensional example of a 



well known limiting case of the n-dimensional Selberg integral [22|, which in general can be 
evaluated as a product of gamma functions. 

To obtain the integral in ([4.4|) we must put z = 1 and p = —1 — (3/2. However we see from 



(A. 4) that L{p,f3) diverges with this choice of p. This is compensated for by a vanishing factor 



in (A^), so the correct procedure is to take the limit p — > —1 — f3/2 (an analogous procedure 
has been necessary in the evaluation of similar integrals occuring in random matrix problems 
|18|). Doing this, and using the gamma function identity 

22^-ir(x)r(x + 1/2) = V^r(2x) 

shows that 

/(-I - P/2, 1, /3) = ^2/5+i(/3/2)2r(-/3/2)r(/3/2 + 1/2). (A.5) 
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